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Abstract. We present a study of the intraday variability be- 
haviour of two samples of x-ray selected BL Lac objects, the 
EMSS and EXOSAT samples consisting of 22 and 1 1 sources, 
respectively. In both samples we were able to detect intraday 
variability in less than 40% of the sources only. The duty cycle 
(the fraction of time, when a BL Lac object is variable) in x-ray 
selected BL Lac objects is 0.4 or less. 

The typical peak-to-peak amplitudes of the variability are 
10%. Typical time-scales and an activity parameter for our vari- 
able BL Lac objects were inferred from structure function and 
autocorrelation function analyses. In only 4 BL Lac objects we 
were able to measure a characteristic time-scale, which was in 
the range between 1.3 and 2.7 days. 

Comparison with our previous study of a complete sample 
of radio-selected BL Lac objects from the 1 Jy catalogue shows 
that x-ray and radio-selected BL Lac objects differ in their duty 
cycle by a factor of 2 and the typical peak-to-peak amplitudes 
by a factor of 3. The observed time-scales are similar. We also 
found that the same mechanism may be responsible for the ob- 
served variability in the x-ray selected and radio- selected BL 
Lac objects. 

The expectations of the various schemes linking x-ray se- 
lected and radio-selected BL Lac objects have been compared 
to our observations. Consistency is found for a scenario, where 
x-ray selected BL Lac objects have on average stronger mag- 
netic fields and are seen under relatively larger viewing angles 
than the radio-selected BL Lac objects. However, the sugges- 
tion that x-ray selected BL Lac objects have decelerating jets 
and radio-selected BL Lac objects accelerating jets can also not 
be ruled out. In any case, any model which explains variability 
on time-scales of days must be able to reproduce the high duty 
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cycle in radio-selected BL Lac objects and the factor of 2 lower 
duty cycle in x-ray selected BL Lac objects. 
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1. Introduction 

BL Lac objects are characterised by strong and rapid variabil- 
ity from radio up to 7-rays, high and variable polarization in 
the radio and optical regime, as well as faint or absent emis- 
sion lines in their spectra. Additionally, in many sources su- 
perluminal motion has been observed. Due to these properties 
it is nowadays believed that their energy output is dominated 
by boosted radiation from a relativistic jet pointing almost di- 
rectly towards the observer. Those objects, whose relativistic 
jets are pointing more randomly onto the sky are thought to be 
the Faranoff-Riley I (FR I) radio galaxies (Urry & Padovani, 
1995 and references therein). 

Whereas the classical BL Lac objects have been identified 
from radio surveys, most of the newly identified BL Lac ob- 
jects have been detected in x-ray surveys. These so-called x- 
ray selected BL Lac objects (XBL) differ from the classical 
radio-selected BL Lac objects (RBL) in many respects. They 
have different q;ro ~ Q^ox spectral indices (Stocke et al., 1985) 
and possibly show different cosmological evolution (Stickel et 
al., 1991, Morris et al., 1991). Additionally, XBL are not as 
strongly variable and not as highly polarized in the optical as 
the RBL (Jannuzi et al., 1993, 1994), they seem to have weaker 
radio cores and extended emission and are less core dominated 
(Perlman & Stocke, 1993). On the other hand, no differences 
have been found among the properties of their host galaxies 
(Wurtz et al., 1996) and their cluster environment at redshifts 
below z = 0.65 (Wurtz et al., 1997). 

Since the detection of "two" classes of BL Lac objects, sev- 
eral suggestions have been made to explain their differences 
and to unify both classes. Ghisellini et al. (1989, 1993) pro- 
posed that XBL and RBL are intrinsically the same objects dif- 
fering mainly by their average viewing angle of their jets to 
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the observer ( < 30° for XBL and < 15° for RBL). Alterna- 
tively, Giommi & Padovani (1994) and Padovani & Giommi 
(1995) suggested that the main difference between XBL and 
RBL is their spectral energy distribution (SED), where XBL 
have their cutoff of the synchrotron spectrum at higher frequen- 
cies (UV/x-ray regime) than the RBL (near-lR/optical regime). 
Since they found transition objects (RBL which had their cut- 
off at higher frequencies and vice versa) they introduced a new 
classification scheme, the so-called LBL (low-energy cutoff 
BL Lacs) and HBL (high-energy cutoff BL Lacs). More re- 
cently, Sambruna et al. (1996) studied the multifrequency spec- 
tral properties of BL Lac objects and explained the differences 
between XBL and RBL by a systematic change of the intrinsic 
physical parameters of the jet such that XBL have higher mag- 
netic fields/electron energies and smaller sizes than RBL. Fi- 
nally, Brinkmaim et al. (1996) proposed that XBL and RBL are 
intrinsically different and either originate from different popu- 
lations or have emission conditions with different physical pa- 
rameters (e.g. accelerating jets in RBL and decelerating jets in 
XBL). In addition to the differences at lower frequencies, RBL 
turned out to be 7— ray emitters at GeV energies, while all of 
the TeV emitting BL Lac sources belong to the XBL subtype. If 
the 7— ray emission is due to inverse Compton scattering, this 
result confirms that the electron population in XBL is domi- 
nated by higher Lorentz factors. 

Intraday variability (IDV ) measurements of BL Lac ob- 
jects in the optical are one tool to probe the physics and geom- 
etry in the innermost part of the AGN (see Wagner & Witzel, 
1995 for a comprehensive summary). Since optical IDV is un- 
affected by extrinsic effects, such as refractive interstellar scat- 
tering, which is important in the radio domain or gravitational 
microlensing, which might only be important on time-scales of 
weeks to months, such studies investigate variations intrinsic to 
the source. With the current availabihty of medium- sized tele- 
scopes in combination with highly efficient CCDs, IDV studies 
of BL Lac objects as faint as 20mag are possible. This allowed 
us to observe large samples of BL Lac objects within a rea- 
sonable time, with an adequate sampling and with errors lower 
than in any other frequency regime. 

Despite the fact that several groups put major efforts in ob- 
serving BL Lac objects on time-scales from hours to months 
(e.g. Fiorucci & Tosti, 1996; Miller & Noble, 1997; Takalo, 
1997) the samples observed are mostly chosen on the basis of 
the brightness of the sources, previously well documented vari- 
ability, or are restricted to specific objects. In most cases XBL 
are missing, since only a few are bright enough to be observed 
with smaller telescopes. 

As a consequence, the optical IDV properties of XBL are 
poorly constrained. Only two XBL - the classical XBL PKS 
2155-304 and Mkn 421 - have been the subject of several de- 
tailed studies, partly during global multifrequency campaigns 
(e.g. Carini & Miller, 1992; Courvoisier et al., 1995; Wagner et 
al., 1997; Pesce et al., 1997; Heidt et al., 1997). These observa- 
tions showed that XBL have a duty cycle close to unity. How- 
ever, in PKS 2155-304 also a quiescent period was observed 
(Heidt et al., 1997). Optical variation on shorter time-scales 



(seconds to hours) has been investigated by Miller & Noble 
(1997) and on longer time-scales by Jannuzi et al. (1993, 1994) 
and Xie etal. (1996). 

In a previous study, we determined the duty cycle, typi- 
cal amplitudes and time-scales of the complete sample of RBL 
from the 1 Jy catalogue (Heidt & Wagner, 1996). The duty cy- 
cle among the RBL is very high (w 0.8), with amplitudes typi- 
cally reaching 30% and time-scales between 0.5 and 3 days. 

Here we present the results of our study of the optical IDV 
properties of two samples of XBL, the EMSS sample extracted 
by Morris et al. (1991) and the EXOSAT sample compiled by 
Giommi et al. (1991). Our goal is to determine the duty cycle, 
typical ampUtudes and time-scales of XBL, to compare our re- 
sults with those derived from our study of RBL and to interpret 
our results in view of the competing models linking RBL and 
XBL. 

This paper is organized as follows: In chapter 2 the samples 
are briefly described, in chapter 3 the observations and the data 
reduction are summarized. Chapter 4 contains the basic results, 
the statistical analysis is presented in chapter 5. In chapter 6 we 
discuss the results of the statistical analysis and compare the 
results with those derived from our study of RBL followed by 
the interpretation. Our conclusions are summarized in chapter 
7. Throughout the paper Hq = 50 km s^^ Mpc^^ and qo = 
is assumed. 

2. The x-ray selected samples 

2.1. The EXOSAT sample 

The EXOSAT sample of BL Lac objects was compiled by 
Giommi et al. (1991) from the EXOSAT High Galactic Lati- 
tude Survey. The 1 1 objects of this sample were selected using 
the following criteria: 

- Equivalent width of emission fines < 5 A (rest frame) 

- 0.2 < aro < 0.55; 0.6 < < 1-5 

- x-ray flux (0.05 - 2.0 keV) > 10"^^ ergs cm'^ s"^ 

- I b I > 20° 

2.2. The EMSS sample 

The EMSS sample of BL Lac objects was extracted by Morris 
et al. (1991) from the Einstein Medium-Sensitivity Survey. It 
contains 22 objects, which were selected according to the fol- 
lowing criteria: 

- Inclusion in the EMSS sample of serendipitous x-ray 
sources 

- Equivalent width of emission fines < 5 A 

- Evidence for a non-thermal continuum in the spectra (Ca II 
break had a contrast of less than 25%) 

- x-ray flux (0.3 - 3.5 keV) > 5 x 10~^^ ergs cm-^ 

- 6> -20° 

Both samples have one object in common (1207.9H-3945). 
Therefore, this object was observed during two runs. Whereas 
for the BL Lac objects of the EMSS sample redshifts are avail- 
able for all sources, only 4 out of 11 BL Lac objects from the 
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EXOSAT sample have published redshifts. The average red- 
shift of the EMSS BL Lac objects is 0.32. When relevant, we 
will set a redshift of z = 0.3 for all BL Lac objects from the 
EXOSAT sample, whose redshift is unknown. 

3. Observations and data reduction 

The observations were carried out during 9 observing runs be- 
tween September 1992 and April 1995 at the Calar Alto Ob- 
servatory in Spain, the ESO in Chile and in Cananea, Mexico. 
The telescopes were equipped with a CCD and an R filter in 
order to perform relative photometry. An overview about the 
observing runs is given in Table L 



Table 1. Observing journal 



Observatory 


Period 


ESO, Danish 1.5m 


September 1992 


Calar Alto, 2.2m 


February 1993 


Calar Alto, 2.2m 


May 1993 


Calar Alto, 2.2m 


March 1994 


Cananea, 2.1m 


April 1994 


ESO, Danish 1.5m 


July 1994 


Calar Alto, 2.2m 


August 1994 


Calar Alto, 2.2m 


December 1994 


Cananea, 2.1m 


April 1995 



In order to create a homogeneous dataset we tried to ob- 
serve each BL Lac object during seven consecutive nights with 
an average sampling rate of 2 hours. Due to various circum- 
stances (weather, technical problems etc.), however, the sam- 
phng and the length of our data trains for the two samples are 
different. 

Each BL Lac object of the EXOSAT sample was observed 
on average 2-3 times per night during 7 nights. Apart from 
EXO 1146H-2455, which could be observed during 3 nights 
only, each object was observed at least during 5 nights. About 
50% of the EXOSAT BL Lacs have been observed at least dur- 
ing 10 nights. 

Our observations of the EMSS BL Lacs were more densely 
sampled. Most of the sources have been observed 4-5 times 
per night for a period of 6 nights, the shortest data train was 
5 nights, with a sampling of once per night (MS 0922.9H-7459 
and MS 1534.2-H0148). About 25% of the sources were ob- 
served during at least 9 nights (c.f. Table 2). 

The integration times varied between 5 and 30 minutes de- 
pending on the brightness of the BL Lac object, the observing 
conditions and the size of the telescope. They were chosen such 
that the signal to noise in the central pixel of the BL Lac ob- 
ject was as high as possible but the count rate approximately 
30% below the saturation limit or the non-linearity limit and 



the count rates of the comparison stars varied no more than a 
factor of two within each campaign. 

The CCD frames were reduced in a standard manner (bias 
subtracted, corrected for dark current, if necessary, pixel-to- 
pixel variations removed using twilight flat-fields). 

In order to carry out relative photometry the count rates of 
the BL Lac object and 5 to 10 comparison stars were measured 
on each frame by simulated aperture photometry. We computed 
the normalized ratio for each pair of objects over the entire 
campaign to construct hghtcurves. By inspecting the Ughtcurve 
of each pair, variable stars were found and rejected from further 
analysis. For the final analysis the Hghtcurves were used which 
included the BL Lac object and the brightest, non-saturated 
comparison star. The errors were estimated from the standard 
deviation (1 a) of the hghtcurve of two comparison stars as 
bright as or fainter than the BL Lac object. They are typically 
in the order of 0.6-2%, in worst cases up to 5% (corresponding 
to 0.06mag). 

4. Variability statistics and amplitudes 

In order to derive the number of BL Lac objects displaying vari- 
ability during the observations we applied a x^-test following 
Penston & Cannon (1970) to each lightcurve. As in our study of 
the radio- selected BL Lac objects from the 1 Jy sample (Heidt 

6 Wagner, 1996) we choose a confidence level of 99.5% as 
cutoff. The x^-test was carried out twice. First we applied this 
test to all hghtcurves. This gave us the fraction of the variable 
and the non-variable BL Lac objects. Since we were interested 
in variability time-scales from several hours to one week, we 
subtracted trends on longer time-scales by fitting a Unear slope 
to each lightcurve of the variable objects. After subtraction of 
these slopes, the x^-test was applied to the residuals. 

6 out of the 11 BL Lac objects from the EXOSAT sample 
displayed variability during the observations, corresponding to 
a fraction of 55%. In only 4 out of the 1 1 sources (36%) we 
were able to detect intraday variabiUty. The fraction of variable 
EMSS BL Lac objects was lower. Here we could detect vari- 
abihty in only 7 out of 22 BL Lac objects (32%), 6 out of this 

7 showed intraday variabiUty (corresponding to 27% for the 
EMSS sample). In Table 2 we list the both samples along with 
their average R band magnitude during the observations, their 
redshift, the number of observations, the maximum time sep- 
aration between two observations and in column 6 the results 
of the x^-tests for a confidence level of 99.5%. A ++ sign de- 
notes the intraday variable, a + sign variable, but on time-scales 
longer as the observing period and a — sign the non-variable 
BL Lac objects. For comparison we give the confidence of the 
X^-test to the residuals after subtraction of a linear slope from 
the Hghtcurves in column 7. The final two columns list the vari- 
ability amplitudes as defined below as well as the measurement 
error and information about the observing period and telescope 
used. 

We checked, whether a lower confidence level, e.g. 90% 
as used by Penston & Cannon (1970), would alter the distribu- 
tion. We ran through the x^ procedure for the lower cutoff Umit 
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as outline above. Whereas we found no change in the classifi- 
cation of the BL Lac objects from the EXOSAT sample, the 
number of variable EMSS BL Lac objects increased. Now, 15 
out of 22 EMSS BL Lac objects (68%) would be classified as 
variable, 8 of them (36%) would display intraday variability. 

Regardless of the cutoff limit used, we conclude that less 
than 40% of all x-ray selected BL Lac objects from our samples 
displayed intraday variability during the observations. 

In order to test, if the variability amplitudes are corre- 
lated with the average brightness of the objects, we calculated 
the variability amplitudes Amp = \/(Ai„ax ^ Ai„i„)2 — 2(7^ 
where Amax and Amin are the maximum and minimum values 
of each hghtcurve and a the measurement errors, respectively. 
They are given for each individual object in Table 2 in colunm 
7. 

The mean ampUtude for the intraday variable EXOSAT 

sources is 9.0±3.8%, for the intraday variable EMSS sources 
10.6±4.5%. Including also the sources variable on longer time- 
scales, the mean amplitude for the variable EXOSAT sources 
is 8.8±3.1% and for the variable EMSS sources 10.0±4.3%. 

In Figure 1 we plot the average brightness of the BL Lac 
objects during the observations against the variabihty amph- 
tudes. For comparison we have also included the "variability 
amplitudes" of the non-variable sources. Obviously, there is no 
correlation between variabihty amphtudes and average bright- 
ness. 
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Fig. 1. Variability amplitudes versus average brightness. The circles 
correspond to the variable, the dots to the non-variable sources. 



5. Temporal characteristics of the variability 

The temporal characteristics of the variability have been stud- 
ied using two statistical methods, structure function analy- 
sis and autocorrelation function analysis. Structure function 
analysis is a powerful tool to measure time-scales from the 
lightcurves, especially when the data samphng is rather inho- 
mogeneous. Additionally, determining the slopes of the struc- 
ture function allows to characterize the underlying physical 
process. The autocorrelation analysis was used to derive typ- 
ical amphtudes per time interval. Since we described the ap- 



plication of both methods for our study of the 1 Jy sample of 
BL Lac objects in detail (Heidt & Wagner, 1996), we will here 
only shortly summarize our adopted strategy. 

First order structure functions (SF) were calculated from 
the lightcurves according to Simonetti et al. (1985). The first 
order structure function is defined as: 



1 ^ 

SF(At) = -^(f(ti + At) 



f(ti))' 



(1) 



The power a of SF oc (At)" characterizes the variability 
and hence the underlying process. If a = 1 shot noise dom- 
inates, towards flatter slopes (a = 0) flicker noise becomes 
more important (Smith et al., 1993). Characteristic time-scales 
show up as maxima in the SF (Wagner et al., 1996). Since we 
were interested mainly on the characteristics of the variability 
on short time-scales and since we wanted to compare our re- 
sults to those derived for the 1 Jy BL Lac objects, we decided 
to study the SF(A t) within the time interval 0.5 days < At < 
5 days. As in Heidt & Wagner (1996) we binned the SF in vari- 
able intervals such that each interval contains the same amount 
of differences AI(ti,t2) with t2 < ti, where the number of 
bins depends on the total number of data points n. 

The autocorrelation function (ACE) of each variable ob- 
ject was computed following the discrete correlation function 
method described by Edelson & Krolik (1988). As in the case 
of the structure function analysis, we had to deal with a rather 
inhomogeneous sampling. Therefore our binsize A t was cho- 
sen to be 0.2 or 0.5 days, depending on the evenness of the mea- 
surements between and during the nights (see Heidt & Wagner, 
1996, section 5.2 for details). 



5.7. Slopes and characteristic time-scales measured from the 
SF 

The slopes a of the SF were measured by fitting a 1^* or- 
der polynomial using the least-squares method to the SF in 
the log(SF)— log(A t) plane. Slopes were determined only in 
those SF, which had at least three bins in the range between 0.5 
and 5 days. This was possible for all 13 variable sources ex- 
cept MS 1757H-7034, where the "best fit" was not acceptable. 
The resulting slopes are given in Table 3. They range from - 
0.6 to 1.6, with a mean value of 0.8 and a dispersion of 0.7. 
The two sources with negative slopes (EXO 1215.3H-3022 and 
MS 1229.2-H6430) are dominated by rapid variability on time- 
scales around one day. Their structure functions are shown in 
Figure 2. 

For two sources of each sample, we were able to measure a 
characteristic time- scale by identifying pronounced maxima in 
the log(A t) - log(SF) plane. These are EXO 1004.0H-3509, 
EXO 1215.3h-.3022 as weU as MS 1221.8H-2452 and MS 
1402.3-H0416. AU sources show similar intrinsic time-scales 
between 1.3 and 2.7 days. 
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Table 2. The x-ray selected samples 



Object 



niR 



dt [days] 



Confidence [%] 



Amplitude [%] 



Date/Telescope 



The EXOSAT sample 



EXO 0556.4-3838 


17.5 




18 


5.20 




99.6 


3.5(0.6) 


09/92 (ESO,1.5) 


EXO 0706.1-^5913 


17.0 


0.125 


13 


9.87 




0.01 


0.1(1.2) 


03/93 (CA,2.2) 


EXO 0811.2-1-2949 


18.5 




11 


9.80 


-1- 


0.01 


9.3(1.2) 


03/93 (CA,2.2) 


EXO 1004.0-1-3509 


19.5 




13 


9.88 




99.9 


9.9(1.5) 


03/93 (CA,2.2) 


EXO 1118.0-F4228 


18.0 




13 


9.90 




79.0 


7.5(2.5) 


03/93 (CA,2.2) 


EXO 1146.9-1-2455 


19.0 




8 


2.12 




5.6 


2.3(1.5) 


03/93 (CA,2.2) 


EXO 1207.9-1-3945 


20.0 


0.615 


9 


10.00 




3.8 


2.2(1.8) 


03/93 (CA,2.2) 


EXO 1215.3-1-3022 


16.0 




19 


5.12 




100.0 


10.2(1.3) 


05/93 (CA,2.2) 


EXO 1218.8-1-3027 


15.0 




29 


9.16 




100.0 


12.3(0.9) 


01/92 (CA,1.2) 


EXO 1415.6-F2557 


17.5 


0.237 


23 


6.16 




34.1 


6.5(1.0) 


05/93 (CA,2.2) 


EXO 1811.7-H3143 


18.0 


0.117 


22 


5.18 


-1- 


3.6 


7.1(1.3) 


05/93 (CA,2.2) 



The EMSS sample 



MS 0122.1-1-0903 


19.5 


0.339 


25 


5.28 




64.9 


6.1(1.7) 


12/94 (CA,2.2) 


MS 0158.5-1-0019 


18.5 


0.299 


15 


4.20 




45.0 


5.0(1.5) 


09/94 (ESO,1.5) 


MS 0205.7-1-3509 


18.6 


0.318 


33 


5.33 




4.5 


4.8(1.5) 


12/94 (CA,2.2) 


MS 0257.9-1-3429 


18.0 


0.247 


35 


5.36 




99.7 


7.4(1.1) 


12/94 (CA,2.2) 


MS 0317.0-1-1834 


17.7 


0.190 


37 


5.36 


-1- 


63.8 


6.8(1.1) 


12/94 (CA,2.2) 


MS 0419.3-1-1943 


19.7 


0.512 


36 


5.38 




87.7 


21.6(4.0) 


12/94 (CA,2.2) 


MS 0607.9-1-7108 


18.2 


0.267 


24 


4.35 




97.0 


7.4(2.5) 


03/94 (CA,2.2) 


MS 0737.9-1-7441 


18.0 


0.315 


24 


4.33 




0.01 


3.0(1.2) 


03/94 (CA,2.2) 


MS 0922.9-1-7459 


19.5 


0.638 


5 


4.08 




93.5 


4.3(1.0) 


03/94 (CA,2.2) 


MS 0950.9-1-4929 


18.6 


0.207 


22 


4.29 




89.7 


4.1(1.0) 


03/94 (CA,2.2) 


MS 1207.9-1-3945 


19.6 


0.615 


7 


5.92 




49.4 


14.1(5.0) 


04/94 (Can,2.1) 


MS 1221.8-H2452 


17.1 


0.218 


23 


4.28 




100.0 


17.5(1.2) 


03/94 (CA,2.2) 


MS 1229.2-1-6430 


16.3 


0.164 


23 


4.29 


++ 


100.0 


6.8(0.6) 


03/94 (CA,2.2) 


MS 1235.4-1-6315 


18.5 


0.297 


17 


6.06 




6.9 


5.6(2.5) 


04/94 (Can,2.1) 


MS 1402.3-1-0416 


16.9 


0.200 


19 


4.23 


++ 


100.0 


14.6(0.8) 


03/94 (CA,2.2) 


MS 1407.9-1-5954 


19.4 


0.495 


24 


9.10 




76.9 


6.7(1.5) 


08/94 (CA,2.2) 


MS 1443.5-1-6349 


19.5 


0.299 


22 


8.07 




88.8 


12.8(2.9) 


08/94 (CA,2.2) 


MS 1458.8-H2249 


16.1 


0.235 


16 


4.19 


++ 


99.5 


6.3(0.7) 


08/94 (CA,2.2) 


MS 1534.2-^0148 


18.3 


0.312 


5 


3.98 




5.4 


0.8(0.8) 


04/94 (Can,2.1) 


MS 1552.1-1-2020 


17.2 


0.222 


25 


7.23 




24.4 


17.9(4.0) 


04/94 (Can,2.1) 


MS 1757.7-1-7034 


18.7 


0.407 


34 


9.16 


++ 


99.9 


9.7(1.5) 


08/94 (CA,2.2) 


MS 2143.4-1-0704 


17.5 


0.237 


33 


9.16 




74.7 


11.5(2.0) 


08/94 (CA,2.2) 



Column 1 gives the name of the sources, column 2 the average brightness (mp ) during the observations and column 3 the redshift. In column 
4 the number of observations are listed, in column 5 the maximum time separation between two observations and in column 6 the result of 
the -tests. A ++ sign denotes the intraday variable, a -i- sign variable, but on time-scales longer than the observing run and a - sign the 
non-variable sources (see text for details). The confidence of the x^-test to the residuals after subtraction of a linear slope from the lightcurves 
is given in column 7. Column 8 contains the variability amplitude as defined in the text. Here the 1 a errors are included in brackets. Column 9 
finally gives the observing dates, the observatory and the telescope used. ESO, 1.5 = ESO Danish 1.5m; CA, 2.2, 1.2 = Calar Alto 2.2m, 1.2m, 
and Can, 2.1 = Cananea 2.1m, Mexico, respectively. 
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Fig. 2. Structure functions of EXO 1215.3+3022 (top) and MS 
1229.2+6430 (bottom). Both sources are dominated by rapid variabil- 
ity on time-scales around one day. This results in negative slopes in 
the range between 0.5 and 5 days. 



Object 


Slope a 


At [days] 


i [%/day] 


EXO 0556.4-3838 


0.59(.12) 




0.32(.49) 


EXO 0811.2+2949 


1.55(.26) 




0.37(.91) 


EXO 1004.0+3509 


1.42(.48) 


2.7(.27) 


0.85(.85) 


EXO 1215.3+3022 


-0.29(.13) 


1.3(.14) 


2.25(.66) 


EXO 1218.8+3027 


0.59(.08) 




0.86(.96) 


EXO 1811.7+3143 


1.44(.22) 




0.36(.57) 


MS 0257.9+3429 


0.55(.14) 




1.14(.47) 


MS 0317.0+1834 


0.92(.19) 






MS 1221.8+2452 


0.38(.15) 


2.4(.25) 


4.94(1.82) 


MS 1229.2+6430 


-0.63(.17) 


2.7(.25) 


2.17(.66) 


MS 1402.3+0416 


1.21(.12) 




2.89(2.04) 


MS 1458.8+2249 


1.44(.ll) 




1.69(.38) 


MS 1757.7+7034 






1.09(.18) 



Table 3. Results of the statistical analyses of the variable BL Lac ob- 
jects. Column 1 gives the name of the sources, column 2 the slopes a 
of the SF, column 3 the time-scales found from the SF and column 4 
the activity parameter I measured from the ACE. 1 a errors are given in 
brackets. The time-scales and activity parameter have been corrected 
for cosmological effects. 



5. LL The activity parameter I from the autocorrelation func- 

tion analysis 

Since we want to measure a typical amplitude per time inter- 
val of a BL Lac object, we calculated the activity parameter 
i = VACF(0)/At from the ACF by fitting a 2""^ order poly- 
nomial to the ACF of each variable object. Therefore the ACF 
was not normalized to unity at A t = 0. ACF(O) and A t were de- 
rived from the maximum and the decorrelation length of the fit 
parabola. With the exception of MS 0317.0+1834 we could de- 
termine I for all variable BL Lac objects. The results are given 
in Table 3 in colunm 4. Except MS 122L8+2452, which had an 
I of « 5 %/day, all variable BL Lac objects displayed a moder- 
ate variabihty of 3%/day or less. 

6. Discussion 

6.1. IDV properties oftheXBL 

Our observation show that only 6 out of 22 (27%) of the EMSS 
XBL and 4 out of 1 1 (36%) EXOSAT XBL displayed IDV. 
Using a less restrict definition of IDV (cutoff limit 90% for the 
test) now 8 out of 22 (36%) EMSS XBL would be classified 
as IDV, the number of EXOSAT XBL displaying IDV would 
not change. Altogether, we conclude that the detection rate is 
less than 40%. Under the assumption that this detection rate 
can be transformed to a duty cycle (the fraction of time, when 
a BL Lac is variable), the duty cycle in XBL would be 0.4 or 
less. 

To our knowledge no systematic study of the IDV be- 
haviour of XBL have been carried out so far. It is hence dif- 
ficult to compare the results obtained with other studies. Xie 
et al. (1996) observed 6 EMSS and 1 EXOSAT XBL dur- 
ing the last years with an irregular sampling, from hours to 
months. All sources have been observed occasionally several 
times per night during two or three subsequent nights. In 5 of 
the 7 sources they found IDV with amplitudes between 0.3 
and O.Smag on time-scales from hours to 2 days. In 3 out 
of these 5 sources (EXO 1218.8+3027, MS 0317+1834, MS 
1402.3+0416) we found also IDV, but no XBL of our samples 
displayed as large variability ampUtudes as found by Xie et 
al. (1996). We found mean amplitudes of w 10 ± 4% in both 
samples, the most extreme case being MS 1221.8+2452 with 
an ampUtude of 17.5%. However, the strong variabihty amph- 
tudes observed by Xie et al. (1996) depend in all cases on one 
single data point in their lightcurves only, so their results must 
be regarded as uncertain. 

Optical photometry and polarimetry on longer time-scales 
(weeks to month) of 37 XBL have been presented by Januzzi 
et al. (1993, 1994). Their sample includes 21 of the EMSS 
XBL (except MS 1443.5+6349) and one EXOSAT XBL (EXO 
1415.6+2559). In spite of the faintness of some sources, which 
made variabihty measurements impossible, they detected vari- 
ability in 12 out of this 22 sources, with amplitudes of always 
less than Imag. The "duty cycle" on time-scales from weeks to 
month of this sample would be « 0.55, which is slightly higher 
as our duty cycle of 0.4. It is interesting to note that their duty 
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cycle of optical polarization (fraction of the time, where a BL 
Lac shows polarization > 4%) of their sample of 37 objects is 
«0.4. 

Optical Variations of 12 XBL on shorter time-scales (sec- 
onds to hours) have been studied by Miller & Noble (1997). 
They found a duty cycle of ~ 0.8 for the XBL they observed, 
which is a factor of two higher as derived from our study 
on longer time-scales. However, they do neither quote their 
sources nor their selection criteria such that a direct compar- 
ison is not possible. 

From our structure function analysis we could derive char- 
acteristic time-scales for 4 XBL only, which corresponds to a 
fraction of 12%. They were in the range between 1.3 and 2.7 
days. The slopes of the SF range from -0.6 up to 1.6 with a 
mean of 0.8±0.7 indicating that in these sources shot noise is 
the most dominating process. From our autocorrelation func- 
tion analysis we could measure a typical amplitude per time 
interval for all variable sources except MS 03 17.0H-1834. In all 
but one source (MS 1221.8h-2452) we derived an i of 3%/day 
which indicates that our XBL show on average modest vari- 
ability only. 

Two prominent XBL - the classical XBL PKS 2155-304 
and Mkn 421 - were subject to several multifrequency cam- 
paigns during the last years (e.g. Carini & Miller, 1992; Cour- 
voisier et al., 1994; Wagner et al., 1997; Pesce et al., 1997; 
Heidt et al., 1997). Whereas the observations generally have 
shown that in both XBL the duty cycle in the optical on time- 
scales of days is close to unity, in PKS 2155-304 also a qui- 
escent period was observed (Heidt et al., 1997). Both objects 
occupy the region typical for XBL in the q;ro — ctox diagram, 
however, they are unusually bright in the radio, optical and x- 
ray domain. Both XBL may be an extreme member of its class. 

We checked for any dependences of the time-scales, ampli- 
tudes or the activity parameter I on intrinsic properties such as 
absolute magnitude or redshift and found none. 

6.2. IDV properties of XBL compared to RBL 

We compare now our results of the IDV properties of XBL to 
those derived in our study of the 1 Jy RBL sample (Heidt & 
Wagner, 1996). 

In Figure 3 we compare the variability statistics of the 
EMSS sample, the EXOSAT sample and the 1 Jy sample. Note 
the clear difference on the number of BL Lacs showing IDV, 
variability on time-scales longer than an observing run (typi- 
cally 7 days) and non- variable sources between XBL and RBL. 
This diagram shows that the rate of occurrence of IDV in both 
classes is different. For our XBL we determined a duty cycle 
of Ri 0.4, whereas a duty cycle of at least 0.8 was found for 
the RBL We checked whether IDV is a coincidence in RBL by 
repeated observations of various subsamples of the 1 Jy sample 
and found always a duty cycle around 0.8 (Wagner et al., 1990; 
Bock, 1994). The duty cycle on time-scales of days differs be- 
tween both classes at least by a factor of 2. 

In order to quantify the difference of the rate of occurrence 
of IDV between RBL and XBL we performed a KS-test to the 



results of the distribution (see section 4). We found a prob- 
ability of^ 2 X 10~^ that the distributions of the have been 
drawn from the same parents. 

Miller & Noble (1997) determined the duty cycle of 12 
XBL and 20 RBL. They found no difference between both 
classes on shorter time-scales (seconds to hours). However, 
they do neither quote their sources nor their selection criteria. 
This is critical for such a comparison. Obviously, the apparent 
discrepancy could be understood, if the samples of Miller & 
Noble consisted mainly of XBL and RBL, which were already 
known to be very active. 



EMSS EXOSAT 1 Jansky 




Fig. 3. Histogram showing the distribution of the sources showing 
IDV (shaded), the sources variable only on time-scales longer than the 
observing run (open) and the non-variable sources (black) among our 
three samples observed. Note the clear difference in the distribution 
between the XBL (EMSS, EXOSAT) and RBL (1 Jansky). 

We also find significant differences in mean amplitudes be- 
tween XBL and RBL. For the intraday variable XBL we de- 
termined a mean amplitude of 10.0±4.1% whereas we found 
mean amplitudes of 28.2±19.7% for the variable RBL. The 
amphtudes in XBL and RBL differ by a factor of three as il- 
lustrated in Figure 4. This confirms that RBL also tend to have 
higher amplitudes on average than XBL in Miller & Nobles 
microvariability study. 

The distribution of the activity parameter I (amplitude per 
time interval) seems to differ among both classes. Whereas the 
variable RBL showed a wide range of I up to 27%/day, the vari- 
able XBL had an i of 3%/day or less (except MS 1221.8H-2452 
with an I of 5%/day). This is shown in Figure 5, where we 
compare the distribution of the I found among the XBL and 
RBL. 

There also seems to be a difference between both classes 
when correlating I with the absolute magnitude of the sources. 
Whereas we found a trend for higher I towards higher Mr in 
RBL, there is no such trend for the XBL as can be seen in 
Figure 6a. However, high luminosity XBL are missing in our 
diagram. In order to test, whether this trend is already present 
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a) 



Fig. 4. Histogram of the amplitudes of the variable XBL (shaded) and 
variable RBL (open). The difference in amplitudes can clearly be seen. 




Fig. 5. Histogram with the distribution of the I for XBL (shaded) and 
RBL (open). 
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Fig. 6. a) Activity parameter I versus absolute magnitude Mr. The 
dots correspond to XBL, the open circles to RBL. The amplitudes 
and time-scales involved in the derivation of I have been corrected 
for cosmological effects, b) Cumulative distributions of the I for XBL 
(solid) and RBL (dashed) having Mr > -27. The KS-test gave a 
probability of 92.5% that the two distributions have been drawn from 
the same parent population. 



in low-luminosity sources, we performed a KS-test on the dis- 
tribution of the i for the XBL and RBL with Mr > -27. The 
cumulative distributions are shown in Figure 6b. The KS-test 
gave a probability of 92.5% that the two distributions are drawn 
from the same parent population. As long as no information on 
high-luminosity XBL is available, we cannot decide, whether 
only RBL show a trend of higher I towards higher Mr or the 
observed trend reflects two intrinsicaUy different types of ob- 
jects. 

In only 4 XBL displaying IDV we could measure a typical 
time-scale of the variations from our structure function analy- 
sis. The range found (1.3—2.7 days) was similar to that derived 
from the 1 Jansky sample of RBL (0.5—3 days). 

The average slopes a of the SF of both classes are similar 
(0.8± 0.7 for XBL and 0.8±0.6 for RBL), implying that the 
same mechanism might be responsible for the observed vari- 
ability in both classes of BL Lac objects. 

Finally, neither class shows a dependence of the amplitudes 
or time-scales on intrinsic properties (e.g. redshift, absolute 
magnitude). 



6.3. Are XBL and RBL intrinsically different? 

Our study of the IDV properties of XBL and RBL has shown 
that both classes differ in many respects. The duty cycle differs 
by a factor of « 2, the variability amplitudes by a factor of w 3. 
Additionally, the typical amplitudes per time interval are mod- 
erate in XBL, whereas a fraction of RBLs show also stronger 
amplitudes per time interval. On the other hand, we found sim- 
ilar time-scales and slopes from the structure function analysis 
among both classes indicating that the same mechanism is re- 
sponsible for the observed variability. In the following we will 
compare these results to the expectations based on the various 
suggestions linking both classes of AGN. 

GhiselHni & Maraschi (1989) and Ghiseflini et al. (1993) 
explained the observed different overall spectral shapes of XBL 
and RBL by a different viewing angle between the relativistic 
jet axes and the observer. Their model is based on the assump- 
tion that the x-ray emission comes from more closely to the 
base of the jet and is more isotropic than the radio emission, 
which originates further away from the baseline and is more 
strongly beamed. They conclude that RBL should be seen un- 
der an inclination angle ^ < 15° and the XBL under an in- 
clination angle 6 < 30°. With this assumption and a Doppler 
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factor 6 scaling as ~ 1/sin^ one would expect higher ampU- 
tudes, shorter timescales and a higher activity parameter I of 
the RBL as compared to the XBL. This is confirmed by our re- 
sults. However, based on their model, one would expect a sim- 
ilar duty cycle among both classes, which we did not found. 
Therefore orientation effects can not account for the observed 
IDV properties alone. 

Alternatively, Giommi & Padovani (1994) and Padovani & 
Giommi (1995) argued that a classification of BL Lac objects 
into XBL and RBL is based rather on the identification in the 
corresponding frequency band as a physical distinction. They 
proposed a new classification scheme, the so-called LBL (low- 
energy cutoff BL Lacs) and HBL (high-energy cutoff BL Lacs), 
where the LBL have the high-energy cutoff of their synchrotron 
spectra somewhere in the IR-optical and the HBL somewhere 
in the UV/x-ray regime. Although most of the XBL are classi- 
fied as HBL and most RBL as LBL there are transition objects 
(XBL, which are classified as LBL and vice versa). In this sce- 
nario simple predictions of IDV characteristics can be made: 
since the spectral slope of the LBL is steeper at optical fre- 
quencies as compared to the HBL one would expect stronger 
amplitudes, similar time-scales and a higher activity parameter 
I for the LBL as compared to the HBL but a similar duty cycle 
for LBL and HBL. 

In order to compare these predictions with our results, we 
sorted our BL Lacs observed in LBL and HBL. All 1 1 EX- 
OSAT and 22 EMSS BL Lacs are HBL and all 34 1 Jy BL 
Lacs except Mkn 501 and PKS 2005-489, which are classified 
as HBL are LBL (Padovani, priv. communication), i.e. we have 
now 35 HBL and 32 LBL. The LBL show indeed stronger am- 
plitudes, similar time-scales and a higher activity parameter I 
as the HBL. However, now 12 out of 35 HBL show IDV, cor- 
responding to a fraction of <40% and 26 out of 32 LBL show 
IDV corresponding to a fraction of 80%. The duty cycle among 
both classes still differ by a factor of 2. Therefore the different 
high-energy cutoff among both classes cannot account for the 
different duty cycle. 

More recently, Sambruna et al. (1996) proposed a system- 
atic change of the intrinsic jet parameters as an explanation for 
the observed different global spectra of XBL and RBL. They 
investigated the multifrequency spectral properties of the com- 
plete EMSS sample, 29 1 Jy BL Lacs and a complete sample of 
flat spectrum radio quasars (FSRQ) from the S5 survey, which 
have all been observed by ROSAT in pointed mode. By com- 
puting composite aox — otx spectra for which positive or neg- 
ative values indicate concave or convex spectra, respectively, 
they found that the EMSS XBL tend to have convex spectra, 
the FSRQ concave spectra and the 1 Jy RBL a mixture of both. 

Based on their results they fitted homogeneous and inho- 
mogeneous (accelerating) jet models to three representative 
(convex, concave and intermediate aox — ctx) multifrequency 
spectra and concluded that a different viewing angle (beaming 
factor) can not account for the observed global spectra alone, 
but requires instead a systematic change of intrinsic jet param- 
eters such as magnetic fields, maximum electron energy or jet 
size. The change should be such that objects having convex 



spectra have stronger magnetic fields, higher electron energies 
and smaller sizes of their jets as objects having concave spectra. 

In order to see, whether there is any trend of the shape of 
the composite spectra with our IDV properties derived, we cor- 
related their ctox — ctx with our amplitudes and I, but found 
no correlation. We also calculated the mean ciox — ctx for our 
variable and non-variable BL Lacs from the EMSS and 1 Jy 
sample and found an aox — ot^ = 0.13±0.61 for the variable 
sources and an aox — Q!x = — 0.24±0.49 for our non-variable 
sources. This shows that the IDV properties are independent of 
aox '^x- 

Nevertheless, their scenario is appealing since it could re- 
produce all of the IDV properties. Under the assumption that 
IDV is most likely originating from perturbances within a rel- 
ativistic jet it might be possible that a varying jet configuration 
such as stronger magnetic fields in the jet of XBL would be 
more resistant against such perturbances. This could explain 
the different duty cycle between XBL and RBL. Additionally, 
a different average viewing angle and hence different beaming 
factors could also explain the different variability amplitudes. 

Different conclusions for the relationship between RBL and 
XBL were drawn by Brinkmann et al. (1996). They favor a sce- 
nario, where either RBL have accelerating and XBL decelerat- 
ing jets or XBL are intrinsically less beamed BL Lac objects. 
Their conclusions were drawn from a multi-frequency study 
of a large number of RBL, XBL and highly polarized quasars 
(HPQ), which have been observed with ROSAT. By comparing 
the x-ray and radio fluxes, the relations between x-ray and ra- 
dio luminosity and by correlating the flux ratios fx-ray/fradio 
and fopticai/fradio of thesc three classes they found indications 
that RBL and XBL are two intrinsically different types of BL 
Lac objects. 

Whereas the scenario that XBL are intrinsically less 
beamed objects can not explain our results (see discussion 
above), the alternative scenario that RBL have accelerating jets 
and XBL have decelerating jets does not provide unique pre- 
dictions on IDV. Additionally, the IDV properties of RBL are 
poorly constrained in x-rays (Heidt et al., in preparation) and 
unknown for XBL in the radio domain. 

Obviously, the main problem for all the models discussed 
above is the different duty cycle between RBL and XBL. One 
possibility to avoid this difficulty could be that most non- 
variable XBL show IDV but below our detection limit. How- 
ever, than there must be a class of BL Lac objects which have 
variability amplitudes a factor of 10-20 lower as compared to a 
typical XBL (0.5-1% versus - 10%). This can hardly be tested 
with the current observing techniques. 

Although several models linking RBL and XBL have been 
suggested in the last years, none of them can fully explain the 
observational facts. This applies also for our observations. Why 
the duty cycle among RBL and XBL differs by a factor of 
two has still not been answered, especially since the physical 
mechanism causing variability is most likely the same in both 
classes. Within the shock in jet model by Marscher & Gear 
(1985) the duty cycle can be transformed in a number of shocks 
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travelling along a relativistic jet. Different duty cycles require 
a mechanism, which injects these shocks with a different injec- 
tion rate. 

7. Conclusions 

We have studied the variability behaviour of the the EMSS and 
EXOSAT samples of XBL on time-scales of days. Our results 

can be summarized as follows: 

- In 4 out of 1 1 EXOSAT sources (36%) and in 6 out of 22 
EMSS sources (27%) we were able to detect variability on 
time-scales of days or less. This impUes that the duty cycle 
in XBL is 0.4 or less. 

- In only 4 BL Lac objects we could measure a typical time- 
scales, which was in the range of 1.3-2.7 days. 

- The typical peak-to-peak amplitudes of the variability were 
~10%. 

- We investigated an activity parameter I from the variable 
BL Lac objects. All but one BL Lac displayed an I of 
3%/day or less. 

Comparison with our study of the complete sample of RBL 

from the 1 Jy catalogue shows that the duty cycle between XBL 
and RBL differs by a factor of two, the typical peak-to-peak 
amplitudes by a factor of three. Whereas the RBL showed a 
wide range of I up to 27%/day, all but one XBL had an I of 
3%/day or less. In both classes we found similar time-scales of 
the variability. We also found similar slopes from our structure 
function analysis for both classes. This implies that the same 
process is responsible for the observed variability. No correla- 
tion of the variability amplitudes or time-scales with intrinsic 
properties have been detected in both classes. 

We compared our results with the various suggestions link- 
ing RBL and XBL. While we can rule out pure inclination ef- 
fects or the different high-energy cutoff of the synchrotron tail 
as being responsible for our observed differences, a combina- 
tion of varying jet configuration (e.g. stronger magnetic fields 
in the jets of XBL as compared to the RBL) and a different 
average viewing angle might explain our observations. A sce- 
nario, where RBL have accelerating and XBL decelerating jets 
can not be tested with our observations. Nevertheless, every 
model explaining variability on time-scales of days must take 
into account the high duty cycle in RBL and the factor of 2 
lower duty cycle in XBL. 
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